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1 Introduction 

The notion of a singularity in general relativity significantly differs from that in other field theories. 
In the absence of a background metric, one has to detect the presence of singularities by showing 
that the space-time is "incomplete" in some sense. In the standard approach to singularities 
(see, e.g., Hawking and Ellis jT^l Ch. 8]), a singularity is regarded as an obstruction to extending 
geodesies. However, this definition does not correspond very closely to ones physical intuition and 
classifies many space-times that have been used to model physically reasonable scenarios as being 
"singular" . Such "weakly singular" space-times have long been used to describe, for example, 
impulsive gravitational waves, shell-crossing singularities and thin cosmic strings. Typically these 
space-times admit a metric that is locally bounded but its differentiability is below C ' (i.e., the 
first derivative locally Lipschitz) — the largest differentiability class where standard differential 
geometric properties, such as existence and uniqueness of geodesies, remain valid. For a recent 
review on the use of metrics of low regularity in general relativity, see [24] . 

This set of problems has stimulated considerations of whether physical objects would be sub- 
jected to unbounded tidal forces on approaching the singularity and was formulated mathemati- 
cally in terms of strong curvature conditions. Unfortunately, it is hard to model the behaviour of 
real physical objects in a strong gravitational field. This led Clarke [2] to suggest that one con- 
sider the behaviour of physical fields (for which one has a precise mathematical description) near 
the singularity instead. According to this philosophy of "generalised hyperbolicity" one should 
regard singularities as obstructions to the Cauchy development of these fields rather than as an 
obstruction to the extension of geodesies. 

However, the weak singularities mentioned above are obstructions if one formulates the Cauchy 
problem for the wave equation in the standard theory of distributions. More precisely, there is no 
generally valid distributional solution concept for the wave equation on a space-time with a non- 
smooth metric. The equation, although linear, involves coefficients of low regularity that cannot 
be multiplied with the distributional solution. 

To resolve this problem in the case of shell-crossing singularities, Clarke [5] introduced a spe- 
cific weak solution concept (called D-global hyperbolicity) to prove unique solvability of the wave 
equation, hence showing that these space-times, indeed, satisfy the conditions of generalised hyper- 
bolicity. On the other hand, Vickers and Wilson [23] used the setting of Colombeau algebras [U[S] 
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to arrive at a valid formulation of the Cauchy problem for the wave equation on conical space-times 
(modelling a thin cosmic string) and showed the existence and uniqueness of solutions in a suitable 
algebra, Q, of generalised functions. Hence they showed that conical space-times are generalised 
hyperbolic or, more precisely, (/-hyperbolic. Vickcrs and Wilson also showed that their unique 
generalised solution corresponds to the "forbidden" distributional solution expected on physical 
grounds (via the concept of association — see Section |2~T1 below) . Their key tool is a refinement 
of the energy estimates for hyperbolic PDEs (see, e.g., [121 Sec. 7.4], [TJ Sec. 4.4]), which makes 
them applicable in the new situation. 

In this paper, we generalise this method to a much wider class of weakly singular space- 
times and prove £/-hyperbolicity for this class. Since our approach is based on regularisation 
of the singular metric by sequences of smooth ones, we must put restrictions on the growth of 
the sequence with respect to the regularisation parameter e. Essentially we shall assume (see 
Section 12. 3|) asymptotic local uniform boundedness with respect to e. Recall that the space- 
times of interest here typically possess a locally bounded metric. In particular, our class includes 
impulsive pp- waves (in the Rosen form), expanding spherical impulsive waves, and conical space- 
times (thereby generalising the results of Vickers and Wilson [25]). 

This work is organised in the following way. In Section [5] we fix our notation, recall some facts 
on the geometric theory of generalised functions and define our class of weakly singular space- 
times (Section I2.3[) . We state our main result in Theorem 13.11 of Section [3) given a point p in 
a weakly singular space-time, there exists a neighbourhood, V, of p such that the initial value 
problem for the wave equation admits a unique solution in G(V). The proof is split into several 
steps: (generalised) higher order energy integrals are introduced and proved to be equivalent to 
suitable Sobolcv norms in Section 3J The energy estimates are provided in Section while some 
auxiliary estimates are proved in Section [6] Finally these results are collected to provide the proof 
of the main theorem in Section [7] We end with some concluding remarks. 

2 Prerequisites 

In this section, we give a precise definition of the class of weakly singular metrics that we are going 
to consider in the sequel. Prior to that, and for the convenience of the reader, we give a brief 
summary of the geometric theory of generalised functions in the sense of Colombeau. Our main 
reference for the latter is [H Sec. 3.2] and we adopt most notational conventions from there. For 
an overview of the use of these constructions in general relativity, we refer to [24| . 

2.1 Geometric theory of generalised functions 

The basic idea of Colombeau's approach to generalised functions [H [S] is regularisation by se- 
quences (nets) of smooth functions and the use of asymptotic estimates in terms of a regu- 
larisation parameter e. Let M be a separable, smooth, orientable, Hausdorff manifold of di- 
mension n, and let X(M) denote the space of smooth vector fields on M. Let (u e )ee(o,il with 
u £ <E C°°(M) for all e. The (special) algebra of generalised functions on M is defined as the 
quotient Q(M) := £m(M)/N(M) of the moderate nets modulo the negligible nets, where the 
respective notions are defined by the following asymptotic estimates: 



£m(M) 



{{u s ) e : VK cc M, Vfc e N 3N e N 
Vi7i, • • • , Vk e X ( M ) : SU P l-^7i • • • ^■n k u e (p)\ 



0(e~ N )}, 



Af(M) 



{(u E ) E ■ VK cc M, Vk, q e N 

VTJi, • • • , Vk e X ( M ) : SU P l-^Ji • ■ ■ U e{P)\ 



Elements of Q(M) are denoted by u = [(u £ ) e ] = (u £ ) e + M{M). With component-wise op- 
erations, Q(M) is a fine sheaf of differential algebras with respect to the Lie derivative with 
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respect to classical vector fields defined by J^ v u := [(jSf^u £ ) e ]. The spaces of moderate resp. 
negligible sequences and hence the algebra itself may be characterised locally, i.e., u G Q(M) iff 
uoij} a g G(il>a(V a )) for all charts (V a ,ip a ), where, on the open set ip a (V a ) C R™, Lie derivatives 
are replaced by partial derivatives in the respective estimates. Smooth functions are embedded 
into Q simply by the "constant" embedding a, i.e., a(f) :— [(f) e ]- On open sets of R™, compactly 
supported distributions are embedded into Q via convolution with a mollificr p G J^(R") with unit 
integral satisfying J p(x)x a dx — for all \a\ > 1; more precisely setting p e (x) = (I / e n ) p(x / e) , we 
have l(w) := [{w * p s )e]- In the case where supp(u>) is non-compact, one uses a sheaf-theoretical 
construction which can be lifted to the manifold using a partition of unity. From the explicit 
formula, it is clear that the embedding commutes with partial differentiation. This embedding, 
however, is not canonical since it depends on the mollifier as well as the partition of unity. A 
canonical embedding of & is provided by the so-called full version of the construction (see [3], 
resp. [lOj for the tensor case). However, since we will model our weakly singular metrics in gen- 
eralised functions from the start (see Sec. 12.21 and the discussion at the end of Section [2TB1 below) 
we have chosen to work in the so-called special setting which is technically more accessible. Note 
that this is in contrast to [25] . 

Inserting p G M into u G Q(M) yields a well-defined element of the ring of constants (also called 
generalised numbers) /C (corresponding to K = R resp. C), defined as the set of moderate nets of 
numbers ((r c ) e G K^ * 1 ! with \r E \ = 0(e- N ) for some N) modulo negligible nets (\r £ \ = 0{e m ) for 
each m). Finally, generalised functions on M are characterised by their generalised point values, 
i.e., by their values on points in M c , the space of equivalence classes of compactly supported nets 
(p e )e £ M^ ' 1 ! with respect to the relation p E ~ p' e :<^> rf/i(p E :Pe) = 0(e m ) for all m, where d^ 
denotes the distance on M induced by any Riemannian metric. 

As is evident from the definitions, all estimates are only required to hold for e small enough, 
that is there exists £o such that for all e < Eq the respective statement holds. However, in order 
not to unnecessarily complicate our formulations we will notationally suppress this fact most of 
the time. 

The C/(M)-modulc of generalised sections in vector bundles — especially the space of gener- 
alised tensor fields Q r s (M) — is defined along the same lines using analogous asymptotic estimates 
with respect to the norm induced by any Riemannian metric on the respective fibers. However, it 
is more convenient to use the following algebraic description of generalised tensor fields 



where 77" (M) denotes the space of smooth tensor fields and the tensor product is taken over the 
module C°°(M). Hence generalised tensor fields arc just given by classical ones with generalised 
coefficient functions. Many concepts of classical tensor analysis carry over to the generalised 
setting [14] . in particular Lie derivatives with respect to both classical and generalised vector 
fields, Lie brackets, exterior algebra, etc. Moreover, generalised tensor fields may also be viewed 
as Q (M )-multilinear maps taking generalised vector and covector fields to generalised functions, 
i.e., as <?(M)-modules we have 



Finally, in light of the Schwartz impossibility result [52], the setting introduced above gives 
a minimal framework within which tensor fields may be subjected to nonlinear operations, while 
maintaining consistency with smooth geometry and allowing an embedding of the distributional 
geometry as developed in [TBI US]- Moreover, the interplay between generalised functions and 
distributions is most conveniently formalised in terms of the notion of association. A generalised 
function u G Q(M) is called associated to zero, u w 0, if one (hence any) representative (u £ ) e 
converges to zero weakly. The equivalence relation u ~ v u — v s» gives rise to a linear 
quotient of Q that extends distributional equality. Moreover, we call a distribution w G 2#'(M) the 
distributional shadow or macroscopic aspect of u and write u ~ w if, for all compactly supported 



Gl{M) = G{M)®X{M), 



(2.1) 



g r s (M) s Lg^ig^MY^liMy-giM)). 



(2.2) 
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n-forms u and one (hence any) representative (u £ ) e , we have 

lim / u € v = (w, u), 



A I 



where {, ) denotes the distributional action. By (|2.ip , the concept of association extends to 
generalised tensor fields in a natural way. 



2.2 Elements of Lorentzian geometry 

A generalised pseudo-Riemannian metric is defined to be a symmetric, generalised (0, 2)-tcnsor 
field g with a representative g £ that is a smooth pseudo-Riemannian metric for each e such that 
the determinant det(g) is invertiblc in the generalised sense. The latter condition is equivalent to 
the following notion called strictly nonzero on compact sets: for any representative (det(g e )) of 
dct(g) we have VK CC M 3m G N : inf p6 if | det(g s )| > e m . This notion captures the intuitive 
idea of a generalised metric as a net of classical metrics approaching a singular limit in the following 
precise sense: g is a generalised metric iff on every relatively compact open subset V of M there 
exists a representative (g £ ) e of g such that, for fixed e, g e is a classical metric and its determinant, 
dct(g), is invertiblc in the generalised sense, i.e., does not go to zero too fast as e — > 0. Note that 
we work exclusively with representatives of generalised metrics that are classical metrics for each 
e. If g is Lorentzian, i.e., there exists a representative which is Lorentzian, we call the pair (M, g) 
a generalised space-time. 

A generalised metric induces a C?(M)-linear isomorphism from Qq(M) to Gi(M). The inverse 
of this isomorphism gives a well-defined clement of Qq(M) (i.e., independent of the representative 
(ge)e)- This is the "inverse metric", which we denote by g" 1 , with representative The 
generalised covariant derivative, as well as the generalised Ricmann-, Ricci- and Einstein tensors, 
of a generalised metric is defined by the usual formulae on the level of representatives. For further 
details see [15]. 

Next, we review the concept of causality in the generalised framework. Let £ £ Qq(M) be a 
generalised vector field on M. Then, by (|2.2p . g(£,£) G Q(M). For functions / G Q{M) we have 
the following notion of strict positivity: 

/ > \/K CC M, 3m G N : inf f £ (p) > e m , as e -> 

and we define time- like, null and space-like for £ by demanding g(£, £) < 0, g(£, £) = 0, respectively 
g(£i £) > 0- (See [IT] for details, as well as for a general account of basic Lorentzian geometry in 
the present setting.) 



2.3 A class of metrics 

We are now ready to define the class of metrics that we will study. Let [M, g) be a generalised 
space-time, and g e a representative of the generalised metric. Let p G M, U a relatively compact 
open neighbourhood of p, and let t : U — > M be a smooth map with the properties that t(p) = 0, 
dt 7^ on U. We assume that there exists an M > with g~ 1 ((it, dt) < —1/Mq, as e — > on U. 
Therefore the level sets of the function t, S T :={<;£[/ : t(q) = r}, are space-like hyper-surfaces 
with respect to the representative metrics, g £ , uniformly as e — > 0. We define the normal covector 
field to these hyper-surfaces er := — dt G 1 (J7) which, via the constant embedding, may also be 
viewed as a generalised covector field on U. We define the corresponding generalised normal vector 
field, £, by its representative £ e G X(U), given, for each e, by cr = g e (£ e , •). We now define the 
generalised function, V, on U by its representative V E : U — > R + , given by 

We will also require the corresponding normalised versions of the generalised normal vector field, 
£ = [(£ e ) e ] = [(£ e /l / e)e]j an d covector field, a = [(3^) e ] = g(£, •)■ Observe that, although cr does 
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not depend on e, the quantities derived from it, i.e., <r^, £ e and £ £ necessarily do, since we are 
dealing with a generalised metric. Using these quantities, one may construct a positive-definite 
metric associated with the generalised space-time (cf. jTTJ Sec. 4]). In particular, we define 

e £ := g e + 25^ <g) a~ e , 

which clearly, for each fixed e, is a Riemannian metric on U . Additionally, the resulting class 
e = [(e e ) e ] defines a generalised Riemannian metric on U ( [TTl Prop. 4.3]). 

We denote by E := Eo the three-dimensional space-like hypersurface through p. Let m be a 
background Riemannian metric on U and denote by || || m the norm induced on the fibers of the 
respective tensor bundle on U . We demand the following conditions. 

(A) For all K compact in U, for all orders of derivative k € Nq and all A;-tuples of vector fields 
Vii - iVk e an d f° r an Y representative (g £ ) £ we have: 

. su VK \\j? Vl ...J? VkEe \\ m = 0(e- k ) ( e -0); 
. su Pif ||^ i ...^g s - 1 || m = 0(e- fe ) (e-0). 

In particular, this implies (for k = 0) that the metrics g e and their inverses g" 1 are locally 
uniformly bounded with respect to s. 

(B) For all K compact in U, we have 

sup||V e a m = 0(l), (e-0), (2.3) 

K 

where V £ denotes the covariant derivative with respect to the Lorentzian metric g e . 

(C) For each representative (g e ) e of the metric g on U, S is a past-compact space- like hyper- 
surface such that 9J+(S) = S. Here J+(£) denotes the topological closure of the future 
emission D+(S) C U of E with respect to g e . Moreover, there exists a nonempty open set 
A <Z and an £o such that 

ac n j+(e). 

£<S 

(Note that we are here following the notation of Friedlander [6]. The set Z)+(E) would be 
denoted /^"(E) in the conventions of, for example, Hawking and Ellis [12].) 

A generalised metric with the above properties will be referred to as a weakly singular metric. 
Some comments arc in order: 

• Condition (A) is independent of the background Riemannian metric, m, chosen on U, and 
may be rephrased in terms of a fixed but arbitrary coordinate system, {x a }, on a neighbour- 
hood of p as follows. With k = 0, Condition (A) states that the components of the metric 
g e and its inverse are locally uniformly bounded on U. In particular, the Lorentzian norm 
of the generalised normal vector field £ e may be assumed to satisfy 

■^J- < V-9e(Ze,e s ) =V S < M . (2.4) 



For k > 0, Condition (A) states that there exists Mk such that 
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• Conditions (A) and (B) imply that the generalised Ricmannian metric, e, obeys the asymp- 
totic condition 

||V s e- 1 || m = (9(l) as e^Q. 

From condition (A), it follows that ||e e || m , 1 1 e^T 1 1| = O(l). (This can most easily be 
deduced from the form of the metric given, below, in (|3.2|) .) Therefore, by the Cauchy- 
Schwarz inequality for the inner product induced by m on the bundle of (2, 1) tensors on 
M, we have 

llVerl = 0(1) as e^O. (2.5) 

ll ^ 1 1 e e 

Similarly, "lowering the index" on £ £ in equation (|2.3[) implies that sup^ ||V 6 er|| m = O(l) 
as e — > 0, where we have again used the fact that ||g e || m = 0(1). Taking the symmetric part 
of V e er implies that ||j^ e g E || = 0(1) as £ — > 0. Finally, again using the Cauchy-Schwarz 
inequality, we deduce that 

11^.^11^=0(1) as e^O. (2.6) 
These estimates will be required in Section [5j 

• Condition (C) is necessary to ensure existence of smooth solutions on the level of represen- 
tatives on a common domain (cf. Step 1 in the proof of Theorem 13.11 in Section [7]) . 

Remark 2.1. Conditions (A) and (B) are given in terms of the e-asymptotics of the generalised 
metric. There is, however, the following close connection to the classical situation. Assume that 
we are given a space-time metric that is locally bounded but not necessarily C < or of Geroch- 
Traschen class [7] (i.e., the largest class that allows a "reasonable" distributional treatment). We 
may then embed this metric into the space of generalised metrics by convolution with a standard 
mollificr (cf. Section [2. ip . From the explicit form of the embedding it is then clear that condition 
(A) holds. 

We recall that in the special version of Colombcau's construction the embedding is non- 
geometric and we could - at the price of technical complications - resort to the full version 
where a geometric embedding is available (as was done in [25]). Nevertheless, in the full con- 
struction generalised functions that arc embeddings of locally bounded functions still display the 
e-asymptotics of condition (A). Moreover our approach using the special version offers more flex- 
ibility: Whenever one succeeds, e.g. by using some physically motivated procedure, to model a 
singular metric by a sequence of classical metrics obeying (A)-(C), then our results apply. 

Condition (B) on the other hand demands somewhat better asymptotics of the derivatives of 
the (0, 0)-component of the metric in adapted coordinates (see also (|3.4[) below). This is a technical 
condition that is satisfied by several relevant examples (see below) . 

As to condition (C), the only part that exceeds the classical condition for existence and unique- 
ness of solutions is the existence of the non-empty open set A. Geometrically, this means that 
the light-cones of the metric g £ do not collapse as e — > 0. In terms of regularisations of classical 
metrics, this condition will always be satisfied if the classical metric is non-degenerate. 

Examples. To begin with we discuss the conical space-times of [33]. They fall into our class 
since estimates (6) and (7) in [25] for the embedded metric imply our condition (A), while (B) is 
immediate from the staticity of the metric. 

The metric of impulsive pp- waves (in "Rosen form") fall into our class. For simplicity we only 
consider plane waves of constant linear polarisation, i.e, 

—dudv + (1 + u+) 2 dx 2 + (1 — u+) 2 dy 2 , 

where u + := uH(u) denotes the kink function. This metric is locally bounded (actually continuous) 
and, since the non-trivial behaviour involves simply the spatial part of the metric, will therefore 
obey Conditions (A) and (B) when embedded with a standard mollifier, or - more generally - if 
we use any other regularisation that converges at least locally uniformly to the original metric. 
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Similarly, in 



metrics for expanding spherical impulsive waves of the form 



2dudv + 2v 2 



dz H — Hdz 

2v 



were studied, where H(z) is the Schwartzian derivative of any arbitrary analytic function h(z). 
Again, this metric is continuous and the non-trivial behaviour occurs in only the spatial directions. 
So we obtain conditions (A) and (B) as for the above case. 

In all of these examples, the discussion at the end of Remark 12.11 imply that Condition (C) is 
also satisfied. 



3 The main result 

We are interested in the initial value problem for the wave equation 

Du = 

u|s = v (3.1) 
3#gu\s = w 

on the subset U of a weakly singular space-time (M, g) (i.e., g subject to the assumptions (A)-(C) 
of Section |2~3")) . Here S := S denotes the level set {q 6 U : t(q) = 0} of the function t :{/—*■ M 
introduced in Section [2. 31 The initial conditions are defined by v and w, which are given functions 
in Note that this, in particular, includes the case of arbitrary distributional initial data. 

We are interested in finding a local solution u G Q on U resp. an open subset thereof. 

A general strategy to solve PDEs in Q is the following. First, solve the equation for fixed e 
in the smooth setting and form the net (u E ) e of smooth solutions. This will be a candidate for a 
solution in Q, but particular care has to be taken to guarantee that the u e share a common domain 
of definition. In the second step, one shows that the solution candidate (u e ) E is a moderate net, 
hence obtaining existence of the solution [(u e ) e ] in Q. Finally, to obtain uniqueness of solutions, 
one has to prove that changing representatives of the data or the metric leads to a solution that 
is still in the class [(it e ) e ]. Note that this amounts to an additional stability of the equation with 
respect to negligible perturbations of the initial data and the metric. 

According to this strategy, given the point p in £ we may, without loss of generality, assume 
that (U, {a; a }), with (x a ) a= o,i,2,3 = {t,x l ) is a coordinate neighbourhood of p, and formulate the 
initial value problem (|3 . 1 [) in terms of representatives on U. To this end, given a representative 
(g e ) e of the metric g, there exist functions hfj,N£ on U such that 

g e = -V e 2 dt 2 + h £ i:j (dx i - N*dt) <E> (dx J - N 3 e dt) . (3.2) 

We further choose representatives (v s ) s , (w s ) s of the data and a negligible net (/ e ) £ on U. We 
then consider the initial value problem 

a s u e = u 

u e (t = 0,x i ) = v e {x l ) (3.3) 
if ? ^ e (t = 0,x l ) = w e {x l ), 

where D £ is the d'Alembertian derived from our particular representative g e , i.e., 

n s u s = \g s \-^d a (\g e \ 1 l 2 g?d b u £ ) 

= -^d 2 t u £ - ^N*dAu e + (h? - ^KN^j didjue - gf r[g £ ] c ab 
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and h l J are the components of the inverse of hfj, g e := detg e , and r[g e ] c a (, denote the Christoffel 



symbols of the metric g e . Note that, by Conditions (A) and (B) of Section |2~B1 the following 
asymptotic estimates hold for the components of the metric in the above coordinate system 



V £ ,h%,Nl = 0{l) 
d a V E = 0(1) 
d a V s ,d a h e ij ,d a Ni 



for all multi-indices a with \a\ > l j 



as e 



(3.4) 



Following the general strategy outlined above, we will prove local unique solvability of (|3.1[) by 
showing that the smooth solutions, (u e ) £ , of (|3.3|) form a moderate net, and hence determine a 
class in Q, and that this class is independent of the choice of representatives of v, w and g. More 
precisely, our main result is the following: 

Theorem 3.1 (Local existence and uniqueness of generalised solutions). Let (M, g) be a gener- 
alised space-time and assume that Conditions (A)-(C) of Section \2.3\ hold. Then, for each p G E, 
there exists an open neighbourhood V on which the initial value problem for the wave equation (|3.1[> 
has a unique solution in Q{V). 

We split the proof in a series of arguments, the core of which are higher order energy estimates. 
To prepare for these, we first introduce suitable energy tensors and energy integrals. 



4 Energy integrals 

By assumption, we have a point p G M and an open neighbourhood of p, Z7, and a map t : U — > R 
with t(p) = such that U is foliated by the level sets of the function t, E r := {q 6 U : t(q) = 
t},t <G [—7,7], for some 7 > 0. Moreover, the level sets E T are space-like with respect to the 
generalised metric g. We now consider solving the forward in time initial value problem for the 
wave equation on U, i.e., with r > (see Figure [1}. 



i 


i t 










P S I 

a 

X 



Figure 1: Local foliation of space-time 

Given p G S = So, let be a neighbourhood of p with the properties that fi C U, and such 
that the boundary of the region fl n {q G U : t(q) > 0} is spacc-lik<u. We denote by S T := S r fl f2 
and by il T the open part of f2 between S and E T . We denote the part of the boundary of fl T with 
< t < t by S n ,r, so that dfl T =S U5 T U Sn, T - 

Notation. In order to simplify calculations, from now on we will adopt abstract index notation 
for (generalised) tensorial objects (sec, e.g., [20 ). In particular, representatives of the metric g e 

lr The existence of such a set, f2, follows from the fact that ||g^~ || m = O(l) as e — » 0. Geometrically, this 
condition means that the collection of timelike directions at a given point is not collapsing to the empty set. 
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and its inverse will be denoted by g e ab and g^ b , respectively, and similarly for the corresponding 
Ricmannian metric e e . We denote the representative of the generalised normal vector field, 
by and the corresponding generalised covector field, <r, by In addition, to simplify the 
notation for tensors we are going to use capital letters to abbreviate tuples of indices, i.e., we will 
write Tj for T^"^J~ with |/| = r, \ J\ = s. Also for I, J of equal length, say r, we write e/j for 
e pi<?i e V2qi ■ ■ ■ e q r pr ■ 

We now use the Riemannian metric e e and the covariant derivative with respect to g e — which 
we have denoted by V e — to define e-dependent Sobolev norms on U. 

Definition 4.1 (Sobolev norms). Let Tj be a smooth tensor field and u a smooth function on U, 
£>0, 0<t<7 and k, j G No. 

1. We define the "pointwise" norm of Tj by 

\\Tj\\l :=e e KL el J T^Tjf 
and the "pointwise norm" of covariant derivatives of u by 

|V«n| 2 := ||V^ • ■ • % 3 u\\l = . . . (v= t . . . V=,u) (v^ . . . V^.n) . 

2. On Q T we define Sobolev norms with respect to VJj resp. partial derivatives by 




where /i e denotes the volume form derived from g e . 
3. The respective "three-dimensional" Sobolev norms arc defined by 




where \x% is the unique three-form induced on S T by /i e such that dt A /j e t = \x e holds on 
S T . Note that although the integration is performed over the three-dimensional manifold S T 
only, derivatives are not confined to directions tangential to S T . 

Observe that, due to the use of a generalised metric, even the norms E depend on 

e. However, due to Condition (A), with k = 0, they are equivalently to an e-independent norm 
derived, for example, from the fixed background metric m. In the following, we will provide 
suitable higher order energy estimates for nets of solutions of the wave equation. These estimates 
are best expressed in terms of energy momentum tensors and energy integrals, which we define 
following [251 Sec. 4]. For the "classical" case, see HH Sec. 7.4], [TJ Sec. 4.4] and [TT] for a recent 
review. 

Definition 4.2 (Energy momentum tensors and energy integrals). Let u € C°°(U) and k € Nq. 
On £1 we define 
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1. the energy momentum tensors by (k > 0) 



2 

2. the energy integrals by 



25* e («) := E / ^'M&M, fc > 0. 

.7=0 •'^ 



(4.1) 



It may be verified, by direct calculation, that the tensor fields T® b ' k (u) satisfy the dominant 
energy condition. Indeed it suffices to observe that, for any future-directed time- like vector field 
U, the expression U a U b — |g e (U, U)g° b defines a Riemannian metric for fixed e. For details, see 
Proposition 3.6 of [TS]. For a generalised formulation of the dominant energy condition, see [T7] , 

Remark 4.3. The energy momentum tensors introduced above are related to the super-energy 
tensors of Senovilla [23] . Omitting indices and e's for the moment, we construct the super-energy 
tensor, S k , of type (0, 2k) (see Definition 3.1 in [23]). Then T k (u) are the (2k - 2)-fold contraction 
of S k with the time-like vector field Theorem 4.1 of [23] then implies that the tensors S k (k > 0) 
satisfy the dominant super-energy property, from which it follows more elegantly that the T k (u) 
satisfy the dominant energy condition. 

Remark 4.4. The energy integrals may be written in the more symmetrical form 



i=a J S T 



where we have defined the volume element jl e T = V e x [f T on S T . In terms of the decomposition of 
the metric given in equation (|3.2[) . p, E T is the volume form on S T defined by the three-dimensional 



metric h e := hfjdx* C3> dx 3 . 



Since the part Sq }T of the boundary of fl is space-like and T" ' J (it) satisfies the dominant 
energy condition, an application of the Stokes theorem yields 

J fl-r J S T ^ So J Sq , t 

J S T J Sq 

where n £ and dS e denote the unit normal and surface element on dfl T , respectively. Hence 
summing over j we have the following energy inequality for each e > and each < r < 7 

E k ^(u) < E k T=0 Ju) + J2 I (6V|Tf •>) + Tf > j (u)Vl£ b ) Me . (4.2) 

j = 

Note that the energy integrals and foliation used here correspond closely to those used in [12l 
Sec. 4.3]. In [25j pp. 1341], due to a different choice of foliation, inequality (|4.2| is replaced with 
an equality. This alternative foliation allows one to work without the explicit use of the dominant 
energy condition, but complicates some of the resulting energy estimates. 

To end this section, we prove the equivalence of the Sobolev norms and the energy integrals. 
Note that this result is the analogue of Lemma 1 in [25] for our class of metrics, and is one of the 
key estimates in our approach. 



5 ENERGY ESTIMATES 



11 



Lemma 4.5 (Energy integrals and Sobolcv norms). 



1. There exist constants A, A' such that for each k > 

A'C\\u\\l !e f < E k ;ju) < AC\\u\\lj 2 (4.3) 

2. For each k >1, there exist positive constants Bk,B' k such that 

( V Hl, £ ) 2 < ^Ei)( a |l-lli,e) 2 (4-4) 

( s im!,j 2 < B k j2-^e\nij 2 (4.5) 

3=1 

Fork = Q we simply have { V \\u\\s T . £ ) 2 = ( a Hls Tie ) 2 - 
Proof. (1): For k = we have 

T^°(u)Ut = -\gfUlu 2 = -^- g£ (Z e ,Z E y = 
hence by l[Z4]). setting A := M /2 and A' := 1/(2M ), we obtain 

< T e ab > (u)£ a g < 

which upon integrating gives the result. 
For the case k > note that 

Hence, we may write 

T^{u)Ul = \v e e?ef<» . . . (V*V^ . . . V^.^XV^ . . . V^u) 

= l -V e \V^u\\ 
Using (A), this implies that 

A'\V^u\ 2 < Tf>i(u)Ut < A\V^u\ 2 , 

which upon summation and integration establishes the claim. 

(2) follows by (A) from the fact that on the compact closure of the metrics e e and S a b are 
equivalent and the Christoffcl symbols and its derivatives are bounded by the respective inverse 
powers of e. □ 

5 Energy estimates 

In this section, we establish the core estimates needed in the proof of our main theorem. 



Proposition 5.1. Let u e be a solution of (|3.3|) on U . Then, for each k > 1, there exist positive 
constants C' k , C k , C' k such that for each < r < 7 we have 



E^ £ {u e ) < Eliu^ + C^UX^ + C'^-^^ f El £ 



(u e )d( 



+Cf / El e {u e )dC. (5.1) 
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Before proving this statement, we draw the essential conclusions from it. Observe that the 
constant in front of the highest order term on the r.h.s. does not depend on e, hence we obtain, 
by an application of Gronwall's lemma. 

Corollary 5.2. Let u e be a solution of (|3 . 3|) on U . Then, for each k > 1, there exist positive 
constants CL, C'^, C^' such that for each < r < 7, 




This statement immediately implies the main result in this section. 

Corollary 5.3. Let u e be a solution of (|3.3j) on U . If, all k > 1, the initial energy (Eq e (u e )) £ is 
a moderate resp. negligible net of real numbers, and (f e ) £ is negligible then 

sup (£* , e (u £ )) E 

0<t<7 

is moderate resp. negligible. 

Proof of \5. 11 Wc begin by estimating the second integrand on the r.h.s. of equation (|4.2|) . Using 
the fact that the energy tensors are symmetric then, by an application of the Cauchy-Schwarz 
inequality to the inner product induced on the tensor bundle Tq(M) by the metric e e , we deduce 
that 

\T?j{u e )VlH b \ < ||Tf <> £ )|| e jV^ b) ||e e = i||T^( We )|| e J|^ E g £ || es . (5.3) 

Equation (|2.6[) implies that there exists a constant K > such that ||jS^ e g e || es < K- In the case 
j = 0, we have 

\\T E ah '°(u e )\\l e = (™) <el h e\ d gTg h e d = u\, 

so 

||Tf>V)lk=^v<% £ | 2 . 

For j > 1, we have 

\\rriab, j /„, Ml 2 „e „s ( 'etc „bd ab cd \ ( a'c' b' d' a'b' c' d'\ 

\\ T e K)IL e =e Qa <e hb , \g e g e - -g s g e I \g e g e - -g e g e I 

x (v c Vjw £ ) (v d Vjit e ) (v c ,Vpu^ (Vd>Vj>u^ el J el' J ' 
= efef (v c V/u e ) (v d V, 7 u e ) (v c >Vjm £ ) (Vd'V./^) e{ J ef J ' 

<A,\V^u e \ 2 , 

where Aj ~ 4 J are combinatorial constants. Letting Ao := 1, we deduce that 

\T^(u e )V^ b \ < l -A 3 K\V^u £ \\ for j > 0. 
Letting A^ :— max_ )= o ! ... ! fc A^, we therefore find that 

< l -A k K(*\\ Ue \\^ £ f (5.4) 




for k > 0. 
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We now consider the first integrand on r.h.s. of (14. 2\i . Beginning with the case k = 1, the 
divergence terms that we require take the form 

VlTf'°(u E ) = -\{Vlg?)ul - Q ff f ) (2u e V> e ) = -u £ V b u £ 

KTf>\u e ) = (g^gf - \gfgf^j (v*v* Ue v^ e + V^ £ V e V^ £ ) 

= V°V»u e V*u e = (D"u e )V b £ u £ = f £ V b £ u £ . 
Inserting this and the k = 1 form of (|5.4|) into (|4.2|) yields 



E\ £ {u £ ) < El £ (u £ )+ [ e^lueUe-u^ + ^KCWu^^) 2 

1 i 

< El £ (u £ ) + Qf (CV^ e )Ve) 2 |/e - U e |Ve) ' + ^1#( V «)* 

1 / 1 1 \ 

< 4 >e (« e ) + Mo QT iv'Viv) 2 I (jf IMVJ 2 + (jf WV=J 2 j 

< El(u £ ) + ^(f (|vWu £ | 2 + | U£ | 2 ) At£ + / ivi^iV+Z IMV) 



2 

< £S llt (« e ) + ^( v ||/.||& r ,.) a +(Mo + iiiiirJ ( v IM f U) z , 

where we have repeatedly used the Cauchy-Schwarz inequality. Now we use (|4.3[) to obtain 

C\K\\^.e) 2 = f T ( V |ke||^. e ) 2 ^C < ^7 f E (J u ^) d C- 
J(=0 A J(=0 

Setting C[ := M /2, C'{ = 0, and C'{' := (M + iii-SQM' = 2M (M + ±iiA") yields the claim 
for fc = 1. 

We now turn to the case k > 1 . We first derive an estimate for 

Z b V e a Tf> k {u e )=h+I 2 + h, 

where we have defined 

h ■= (g?& ^g c A (Ke 1 /) (V*VK)(V £ d V> £ ) 
h ■= -2ei J mv*ju s ) (GflT'Vf^Vf^ 

h ■■= el J (^V> £ ) ( 5 rV*V*VK) 

The strategy is, again, to remove the terms involving derivatives of order k + 1 using the wave 
equation. This requires interchanging the order of covariant derivatives, and therefore introduces 
additional curvature terms. We now calculate the moduli of the terms Ii, I2, I3 separately. 
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We begin by estimating |7i|: 


\h\ = 




< 


(gTtt 







&9e 



cd 



V £ e /J 



V £ e /J 



• ||(V^ e )(V e d V> £ 

•|vw ue | 2 , 



where the inequality in the second line results from applying the Cauchy-Schwarz inequality to 
the inner product induced on the tensor bundle T^ k (M) by the metric e e . The square of the first 
term may then be evaluated as 



1 



9eTs ~ ^e9f ) K^/ 



1 



a r cd 



a c fd ! ^_fa! c d' 



fa" ) (V a ,ef J ' 

/' (y%4 J ) {^i4 J ') = \%\\ 

We now note that, by Condition (A) and equation (|2.5|) of Section HOI we have, 



= 116111 e a /eh,< 



Il'tJJ' 



Iv^ll 2 



|V^ J |I =0(1), (e-0). 



In particular, on each compact set there exists a positive constant, C k , such that ||V|e^ J || e < C k , 
as e — > 0. Therefore, we have the following estimate for I\: 



\h\ < C k ■ ||&|l ■ |V( fc V| 2 < C fc A/ • |VW« 6 



(5.5) 



locally, as e — > 0. 

Next we turn to /2- We then have 



2ei J (VSVX) (Cff £ cd Vf a V^Vfu £ 



. J./ 



(V|V> E )(Cef tV=,VaVf« e )| 



[V*,V*]V?u £ 



where the equality in the first line follows from skew-symmetry in a, c, and the inequality on the 
third line follows from the Cauchy-Schwarz inequality Moreover, from Condition (A), we have 
the following estimates for the curvature on compact sets 



e r>e c 



V e R 



< 



I > 0. 



(5.6) 



Using this estimate with / = and the Ricci identity, we deduce the existence of a combinatorial 
constant N k depending only on k such that 



<N k ^\V^u e \ 



Hence, we have 



\h\ < N k ^ ll&IL • |V«u.| • |V( fc -^ £ | < (|V^ £ | 2 + ±\^u £ A . (5.7) 
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on compact sets. 

Finally, we consider the term 73. We then have, by the Cauchy-Schwarz inequality 
\h\ = \e{ J (4 d VlV> £ ) ( 5 rV^V} U£ )| < ||^V £ d V^ £ || e£ • ||. 9 rV e Q V^V^ £ || ee 

<P fc ||ee|| ee -|Vi fe ) Ue |.||5rV^Vf Ue || ee , 

where P k is a combinatorial constant depending only on k. Again using the Ricci identities, and 
the fact that u e is a solution of (|3.3p . we may write 

fc-i 



where 7^.£ fc llJ '^u e denotes a linear combination of contractions of the (k — j — l)'th covariant 
derivative of the Riemann tensor with the j'th covariant derivative of u £ . A second appeal to (|5.6| 
implies that on each compact set there exists a constant Gk such that 



< 



G k 

£ k-j+ 



r |V«u e |, (e-0). 



Wc therefore have 

|/ 3 | < p k ii^iu ■ iv^vi [ \^t^h\ + G k x ^-|V 



fc-i 



< 



|v^/ e | a + ^E^™|v 



Putting together ()5.5|) . (|5.7|) . and (|5.8|) . we have 



W u I 2 



— 1 'V^u ' 2 



(5.8) 



for positive constants afc,/3fc,7fc. Summation over = l...m and integration yields positive 
constants 5 m , /3 m , 7 m such that 



m ~ 



m— 1 1 

< 5 m ( v ||^||^, e ) 2 +^( v ||MII?-i) 2 +7 ro E F 2(i +m - j) ( V 



I^IIL.e) 2 



i=i 



On substituting this inequality and (|5.4[) into equation (|4.2p . we deduce that 

7^ e (u E ) < E™ £ {u e ) (5.9) 
1 \ ~ m_1 1 

( V |HII? T , E ) 2 + /U V IIMI0 2 +^E ^ci^ T ( v IKlli,e) 2 - 



+ a 



As in the case with k = 1 , we may use Lemma 14.51 to write 

( v iHii >e ) 2 = r c\\u E r S( . £ ) 2 d(<±j f T Ei e (u e )d(, 

for j = 1, . . . , m. Substituting these relations into (|5.9|) yields the inequality (|5.1j) . with C' m := (3 m , 
C£ := y m /A' and := (5 m + ±i m X)M'- □ 

Remark 5.4. As can be seen from the expression for VqT" 6,0 (u 6 ), there is no estimate of the 
form (|5.ip for E® e {u e ). However, E® e (u e ) is estimated in terms of E^ £ (u e ), with fc > 1; a fact 
that is implicit in Proposition [57TJ 
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6 Auxiliary estimates 

In this section, we complement the energy inequalities derived in Section [5] with estimates that 
allow us to utilise the former in the proof of the main result. In particular, we shall prove that 

(i) suitable bounds on the initial data give suitable bounds on the initial energies Eq £ (u £ ); 

(ii) suitable bounds on the energies E\ £ (u £ ) give suitable bounds on the solution u £ . 

The existence as well as the uniqueness part of the proof of the main theorem will then use (i) 
combined with Corollary [O] and (ii) to establish moderateness resp. negligibility of the candidate 
solution. 

Lemma 6.1 (Bounds on initial energies from initial data). Let u £ be a solution of (|3.3p . If 

(v £ ) £ , (w e ) s are moderate resp. negligible, then the initial energies (Eq e (u £ )) £) for each k > 0, are 
moderate resp. negligible nets of real numbers. 

Proof. The estimates for the spatial derivatives d x ii . . . d x i k u £ (0, x 1 ) = d x ii . . . d x i k v £ (x l ) are im- 
mediate. To estimate dtd x ii . . . d x % k u s (Q, x % ), we rewrite the initial conditions in equation (|3.3p in 
the form 

u s (t = 0,x l ) = v^) 
d t u e (t = 0,x l ) = w £ (x l ), 

where we define w £ := V £ w £ — N £ d x iV £ . It is straightforward to show, using the asymptotic 
estimates (|3.4p . that (v £ ,w £ ) being moderate resp. negligible implies moderateness resp. negligi- 
bility of (v £ ,w £ ). Therefore moderateness resp. negligibility of (v £ ,w £ ) implies moderateness resp. 
negligibility of d t d^ t . . . d x i k u £ (0, x l ) = c^i ■ ■ ■ d x i k w e (x l ). 

The estimates for higher (mixed) time derivatives follow inductively by rewriting the wave 
equation in the form 

dhe = -V £ 2 U + ^N%diu s - (h? - ^N* e N^ did jU£ + gf T[g s ] c ab 

and using again the estimates (|3.4p for V £ , N £ , h l J as well as f £ , v s , w s . □ 

Lemma 6.2 (Bounds on solutions from bounds on energies). Form > 3/2 an integer, there exists 
a constant K and number N such that for all u £ C°°(fl T ) and for all £ £ [0, t] we have 

sup \d x a, 1 ---d x a l u{x)\<Ke- N sup E" l + l (u). 
xen T o<c<r 

Remark 6.3. Note that the statement is for all u £ C°°(il T ). In the proof of the main theorem, 
we will apply it to a solution, u £ , of the wave equation. 

Proof of \6. 2\ First we combine the standard Sobolev embedding theorem on S T with the fact that 
by assumption (A) the metric and hence the volume is 0(1) to obtain for m > 3/2 

sup | u (a:)| <K d \\u\\^ £ . (6.1) 

Then we successively apply f|4.5[) and (|4.3p to obtain 

sup |«(a;)| < e- n E™ £ (u). 

Taking the supremum over £ G [0, r] on the right hand side gives the result for I — 0. To prove 
the general result, we replace u by the respective derivatives. In some more detail, note that time 
derivatives are not covered by the Sobolev embedding theorem since they are transversal to S T , 
i.e., we have to replace (|6.ip by the estimate 

sup \d Pl . . . d Pk d s t u\ < K d \\dtu\\™+ k £ < K d \\u\\%+* +s , 

xeSc 1 1 

where the last inequality holds because the norm d \\ ||™ E , in addition, contains time derivatives. 

□ 
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7 Proof of the main theorem 

We finally prove the main result by putting together the estimates achieved so far. 
Proof o/EH 



Step 1: Existence of classical solutions. Due to assumption (C), classical theory provided us with 
smooth solutions for fixed e. More precisely, by [51 Theorem 5.3.2], for e fixed there exists a unique 
smooth function u £ solving (|3 . 3[) on A C n e < £0 Jt(^)- Without loss of generality, we may assume 
that f2 7 C A. 

Step 2: Existence of Q -solutions (moderateness estimates). We show that the net (u £ ) £ of Step 1 
is moderate on f2 7 : Moderate data (v £ ) £ , (w e ) e translate, by means of Lemma [6.11 to moderate 
initial energies (Eq e (u s )) £ for each k > 1. Moreover, by means of Corollary [531 moderate initial 
energies translate to moderate energies (E* E {u e )) £ {k > 1) for all < t < 7. Finally, it follows 
from Lemma [6.21 that moderate energies (E^ £ (u £ )) £ (k > 1, < r < 7) imply moderateness of 
(u £ ) £ . Hence u := [(u £ ) £ ] is a generalised solution on fl T of the i.v.p. (|3.3p . 

Step 3: Uniqueness of Q-solutions (negligibility estimates). We are left with showing that the 
solution u does not depend on the choice of representatives of {f £ ) £ , {v £ ) £1 (w £ ) £ , and (g e ) e of 
/ = 0, v, w, and g. Leaving the latter for Step 4, we observe that, to show independence of the 
choice of representatives of /, v, and w, it suffices to prove that if (v £ ) £ and {w £ ) £ are negligible, 
then the corresponding solution (u £ ) £ is also negligible. To establish this claim we argue as in 
Step 2 but using the negligibility parts of Lemma 16.11 and Corollary 15.31 We then observe that 
negligibility of the energies in Lemma 16.21 implies negligibility of (u £ ) £ . 

Step 4-' Independence of the representative of the metric. Wc finally prove independence of the so- 
lution on representatives (g £ ) e of the metric. So let (g e ) e be another representative of g. Denoting 
the corresponding d'Alcmbcrtian by D £ wc consider the initial value problem 

D £ u £ = f £ , 

u £ (t = 0,x l ) = v £ {x l ), (7.1) 
d t u £ {t = Q,x l ) = w £ {x l ). 

By Step 2, there exists a moderate net of solutions (u £ ) £ of (|7.1[) . and we only have to show that 
its difference with the unperturbed solution, (u £ ) £ :— (u £ ) £ — (u £ ) £ , is negligible on Q T . This 
difference is a solution of the i.v.p. 

&u £ = f £ -U e u £ 
u £ (t = 0,x i ) = (7.2) 
d t u £ (t = 0,x i ) = 0. 

In view of Step 3, we only have to show that f £ — \3 E u £ is negligible. To this end, we write 

f £ - &u £ = (f £ - D £ u £ ) + (D £ u £ - D £ u £ ) = D £ u £ - D £ u £ , 

where we have used the fact that {u £ ) £ solves (|3.3p . Therefore, the problem is reduced to showing 
that (d £ u £ — D £ u £ ) £ is negligible. This, however, is clear since □ is a well-defined differential 
operator on Q. □ 



8 Conclusion 

We have proved unique local solvability of the wave equation for a large class of metrics of low 
regularity in the Colombcau algebra of generalised functions, hence establishing t/-hyperbolicity of 
these space-times in the sense of Vickers and Wilson [25] . (This, in itself, is a slight modification 
of Clarke's notion of generalised hyperbolicity [3j, in the sense that we now consider solvability 
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in Q rather than 3)'.) The essential assumption on this class of metrics is local boundcdness: 
in particular, it includes conical space-times, and therefore generalises the results of Vickers and 
Wilson |25j . Our class of metrics also includes non-static examples such as impulsive pp- waves 
and expanding spherical impulsive waves. 

Finally, we remark that the regularity assumptions (A) and (B) on the metric may be relaxed 
slightly. Indeed, we can replace the 0(l)-asymptotics for the zeroth order derivative of the metric 
in Condition (A) as well as in Condition (B) by the condition that these quantities be 0(log(l/e)). 
(This corresponds to generalised Holder- Zygmund regularity of order zero of the metric as defined 
in [13] •) Under these conditions, the constants A and A'^ 1 in (|4. 3|) of Lemma 14. 51 as well as of C(" 
in Proposition 15.11 and Corollary 15. 21 have a growth behaviour of C>(log(l/e)). However, Corollary 
15.31 remains unchanged since the 0(log(l/e))-growth together with Gronwall's lemma still yield 
moderateness resp. negligibility estimates. Therefore, given a classical metric which we regularise 
(cither by convolution with a mollifier or by any physically motivated procedure) subject to these 
weaker asymptotic conditions, then our existence and uniqueness result still holds. 

From these considerations, we also see that it is hard to imagine how the regularity assumptions 
for the metric could be further relaxed within our framework. 
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